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SUMMARY 


The objective of this paper is to develop a more accurate procedure for 
the determination of the inelastic behavior of structural components. For 
this purpose, the actual stress-strain curve for the material of the structure 
is utilized to generate the force-deformation relationships for the structural 
elements, rather than using simplified models such as elastic-plastic, 
bilinear and trilinear approximations. 

Force-deformation curves in the form of universal generalized stress- 
strain relationships are generated for beam elements with various types of 
cross sections. In the generation of these curves, stress or load reversals, 
kinematic hardening and hysteretic behavior are taken into account. 
Intersections between loading and unloading branches are determined through 
an iterative process. 

Using the inelastic properties determined in this study, the plastic 
static response of some simple structural systems composed of beam elements 
is computed. Results are compared with known solutions, indicating a consi- 
derable improvement over response predictions obtained by means of simplified 
approximations used in previous investigations. The application of this 
procedure to the dynamic load case is currently in progress. 


INTRODUCTION 


Structural systems analyzed and designed for traditional loads and 
materials have been observed to undergo inelastic deformations when excessive 
load conditions are experienced. It is, therefore, an established fact now 
that inelastic deformations do occur in structures and are considered in the 
analysis in order to produce more economical and safe designs. For example, 
a generally accepted philosophy in the seismic analysis and design of struc- 
tures is that *a structure should remain elastic during earthquakes of small 
intensity that occur frequently; it should undergo limited plastic deformations 
during earthquakes of moderate intensity; however, it may undergo large plastic 
deformations but without major collapse during earthquakes of relatively high 
intensity that occur infrequently. 

*This study has been partially supported by the National Science Founda- 
tion Research Grant No. PFR-79-16263 . 
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At the same time, advances in naval, aerospace and nuclear reactor 
technology has led to the use of new materials such as stainless steel, alloys 
of aluminum and nickel, reinforced plastics, etc. The stress-strain curves 
for these materials are generally nonlinear. Therefore, an economical design 
of structures composed of such materials requires an accurate knowledge of the 
magnitude and distribution of the stresses and strains, as well as, the 
displacements. In all cases the effect of nonlinearities must be considered 
in the analysis. 

In the study of the inelastic behavior of structures various idealizations 
to the actual stress-strain curves or force-deformation relations have been 
employed. The most extensively used model is the elastic perfectly plastic 
representation, principally due to its simplicity. When unloading occurs, 
this model neglects the strain hardening and Bauschinger effects. In general 
it produces conservative results and is mostly suitable for mild steel 
structures. Bilinear models with a nonzero slope for the inelastic branch 
have also been used widely. These models allow for the consideration of strain 
hardening effects (both isotropic and kinematic) due to loading and unloading 
cycles arising from static and dynamic loads. A trilinear model has also been 
used to simulate the stress-strain relationship of the material under static 
loads. 


Simplified models often perform satisfactorily in predicting the inelastic 
behavior of special classes of structures. However, for general types of 
structures, an accurate representation of the stress-strain or force- 
deformation properties are needed for both the loading and unloading branches 
in the form of curvilinear or multilinear (piecewise linear) relations that 
follow as closely as possible the actual behavior of the system. 

In the present study an accurate procedure is considered for the determi- 
nation of the inelastic behavior of structural components. For this purpose, 
the actual stress-strain curve for the material of the structure is utilized 
to generate the force-deformation relationships for the structural elements, 
rather than using simplified models, such as elastic-plastic, bilinear and 
trilinear approximations. Applying this process to frame type structures, 
force-deformation curves in the form of universal generalized stress-strain 
relationships are generated for beam elements with various types of cross 
sections. In the generation of these curves stress or load reversals, more 
realistic strain hardening properties and hysteretic behavior are taken into 
account. Intersections between loading and unloading branches are determined 
through an iterative process. Based on the rather accurate force-deformation 
relationships of the individual elements, the governing equations for the 
structural system are established and used to compute the inelastic response 
of the structure. 
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GOVERNING EQUATIONS 


Application of the Hamiltonian Principle to discrete systems in the 
context of the finite element method yields the basic dynamic equations 
governing the behavior of the structural systems. In matrix form these 
equations can be expressed as (ref. 1) 

[m] {q} + ( [K] + [K G ]){q} = {F} + {F°} (1) 

where [m] = consistent mass matrix of the structural system. 

[K] = elastic stiffness matrix of the structural system. 

[K ] = geometric stiffness matrix of the structural system. 

{q} = vector of displacements at the structural degrees of freedom. 

{q} = vector of accelerations at the structural degrees of freedom. 

{F} = generalized nodal force vector corresponding to externally 
applied loads. 

{F°} = equivalent generalized nodal force vector due to plastic 

strains, computed in accordance with the initial stiffness 
method. 

In case of static loading the above equations take the form of 

C [K] + [K g ]) {q} = {F} + {F°} (2) 

PLASTICITY RELATIONS 


Stress-Strain Curve 

An experimentally determined virgin curve of the material is in general 
curvilinear. Starting with the experimental stress-strain data analytical 
expressions can be obtained to represent this data. Such expressions can be in 
the form of algebraic or other types of polynomials, exponential functions, or 
the widely used curvilinear relationship known as the Ramberg-Osgood approxima- 
tion represented by 



(3) 
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in which e and o are the unit strain and unit stress, respectively, E repre- 
sents the modulus of elasticity, a y denotes the yield stress of the material 
and D and R are real constants to be determined. However, this relationship 
is not explicit in stresses and, therefore, numerical procedures are needed to 
find the stresses corresponding to given strains. On the other hand, if the 
stress-strain data cannot be represented by an analytical expression, the 
curve is approximated by a series of line segments given by 


cr = 


CT ± + k ±+i (e - G i ) 


e . < e < e . , . 
x — x+1 


(4) 


where, (a. ,e.) are the stress-strain values at the beginning of the ith segment 
and k^+2 is the slope of the line segment between the points (a.,e.) and 
(a i+ ^,e i+ ^) expressed as 


k 


i+1 



, i = 0,1,2, 


(n-1) 


where n is total number of segments used to approximate the curve. 


(5) 


Moment-Curvature (Force-Deformation) Relationship 

Figure 1 shows a general elasto-plastic doubly symmetric cross section of 
a beam member. A linear strain distribution over the depth of the cross 
section up to ultimate behavior is assumed. Tensile stresses are considered 
positive and curvatures causing positive strains at bottom fibers are also 
positive. The x- and y-axes are the principal axes and the z-axis is the 
geometric centroidal axis of the cross section. The bending moment acting on 
the cross section is the sum of moments of the stresses acting on the cross 
section about the geometric centroidal axis, i.e. 

M = J o ydA (6) 

A 

in which dA is an element of area. The integration is carried over the elas- 
tic and plastic parts of the area. Usually the stress-strain curve of the 
material cannot be represented explicitly for stresses; therefore, a numerical 
integration procedure is employed. For this purpose the section is divided 
into a series of rectangular slices, and the contribution of each slice to the 
moment acting over the cross section is found. Before indicating the details 
necessary for accomplishing this, certain definitions need to be established as 
follows. If £y, e u and a y are the yield strain, ultimate strain and yield 
stress, respectively, as observed in a tension test on a material and h is the 
distance to extreme fibers from z-axis, then 
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M = a I /h = a S 
y y z y 

where My, I z , S are the yield moment, moment of inertia and section modulus 
of the section, respectively. The yield curvature $ is obtained from 


( 7 ) 


= M /E I 
y z 


( 8 ) 


Based on the assumption of a linear strain distribution over the cross section 
of the member, the ultimate curvature is given by 


= e /h 
u u 


(9) 


In the present study the curvature range between <py and <j> u is divided into 
a suitable number of intervals to give enough data points to fit a curve. 

The strain e. at the center of ith rectangular slice, assumed to be uniformly 
distributed over the slice, corresponding to a curvature value <f) . is given by 


e . 

i 


= <J> . y . 
J J x 


( 10 ) 


where y^ is the distance to the centroid of the ith rectangular slice from the 
centroidal axis of the cross section. The strain in the ith slice due to 
curvature cj> j , at the j th discrete point along the member, is used to determine 
the stress Oj_ assumed to act uniformly over the entire slice, from the stress- 
strain relationship as discussed previously. If is the area of the ith 
slice, the force acting on the ith slice is given by 


F. = a. A. 

1 IX 

The moment of the force about the neutral axis is obtained from 


(ID 


M. = F.y. (12) 

l i y l 

By summing over the total number of slices used to model the cross section, the 
total internal moment Mj on the cross section corresponding to a curvature ( p . 
is expressed as 


NS NS 

M. = I M. = Z a . A. y . 

J i=i 1 ±=i 1 1 1 


(13) 


where NS is number of slices and j represents a typical discrete point along 
the member. Similarly the total internal axial force P at the cross section 
due to a prescribed strain field e is obtained from 


NS 

P = Z a. A. (14) 

±=i 1 1 
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Equations (13) and (14) directly yield force P and moment M corresponding to 
a prescribed deformation field. However, to get the deformation corresponding 
to a given load history^ equations (13) and (14) have to be solved iteratively. 
The moment curvature data generated through the use of equation (13) is then 
presented explicitly for curvature by a polynomial of the type 

N 

<f> = M/EI + Z C. M 1 (15) 

1=1 1 

where C. are the polynomial coefficients to be determined and N is degree of 
the polynomial. Alternatively, an equation of the Ramberg-Osgood type 

<}> = £ 1 + D(M/My) J (16) 

can be fitted in which D and R are the real constants to be determined. 


Normalized Moment Curvature Relationship 

The moment-curvature data generated from Eq. (13) is normalized by 
using the quantities M^ and (p^ 9 so that 

M = M/M (17) 

y 

* = <P/<t> y (18) 

where M and <f> are normalized moment and normalized curvature, respectively. 

The normalized curvature is then separated into its elastic and plastic 

components (p and cj)p, respectively, so that solving for <f> 
e p 

+ = * - * e (19) 

Since (p = M/EI, based on equation (17) the normalized elastic curvature is 
written in the form 

<j> = M 

Y e 

Based on the normalized plastic components of curvature computed from 
equation (19), a polynomial of the form 
N 

<f> = E d. M 1 

P 1=0 1 

is then fitted to approximate the data, where the d.'s are constants to be 
determined through regression analysis. Alternatively, an equation of the type 


( 20 ) 


( 21 ) 
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<b = D M R 
P 


( 22 ) 


is also used to represent the data. The real constants D and R are determined 
by taking the logarithmic form of Eq. (22), so that 

In i * In D + R In M (23) 

P 

This is an equation of a straight line. Constants D and R are determined from 
linear regression analysis using the data for In <j> and In M. 

Using equations, (19) , (’21) , and (22) , the following moment-curvature 
relations can now be formulated 

cj> = M , M < 1 (24) 


N 

<j> = M + E d. M 1 , M > 1 (25) 

1 

i=0 


<()=M + AM R = M [1 + AM 11 1 ] , M > 1 (26) 

Slope of Generalized Stress-Strain Curve 


A universal stress-strain curve usually represents the relationship 

— P — 

between <j) and M expressed by Eqs. (21) or (22). Rewriting these equations in 
differential form and solving for the slope K 2 yields 


T 7 _ dM 

K 2 -P 
dcj> 


N 

£ 

i=l 


(27) 


i d. M 1) 
1 


and 


K 2 = 


dM 

d^P 


AR M (R - 1} 


(28) 


-P 


If a smooth analytical expression to fit the <j> - M data is not possible with- 

in tolerable limits of accuracy, the slope of the universal stress-strain 
curve is obtained as a series of approximate tangents drawn at discrete points 
representing the data, i.e. 
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(29) 


(K 0 ), = 


M - M. 
l+l 1 


2 i TP TP. 

♦ i+1 " + 1 


where i and i+1 are typical discrete points. 


In eqs. (21), (22) and (24)-(29) the cross section is assume^ to be 
subjected to a bending moment about the centroidal axis. However, in practical 
cases, structural elements are generally subjected to stress resultants acting 
in different directions. To extend the applicability of the above procedure, 
it is assumed that similar relationships exist between the effective stress 
f* and effective strain 0* in a multi-axial stress case. The function f* in 
normalized form is also identified with the yield function. 


Yield Function 


When stress resultants in normalized form (instead of unit stresses) are 
considered to be generalized stresses in the context of plasticity theory, 
yielding at any section of a member is then assumed to occur when the critical 
combination of generalized stresses initiate inelastic deformations at that 
section. The yield function is expressed by an equation of the form 

f*(Q_.) = Y (30) 

where Q. are the generalized stresses and Y is the initial yield value. To 
make the yield function independent of the cross section, the yield function 
equation is derived in terms of normalized (dimens ionless) force parameters. 

A force component is normalized by its corresponding characteristic value 
(usually the value at first yield). In normalized form equation (30) takes 
the form of 


f*(Q ± ) = 1 (31) 

For a space frame member the cross section is subjected to a generalized force 
vector S having 6 components > so that 


{S} = 



M M , M } 


x 


where 

P = axial force 
x 

V , V = direct or transverse shear forces 

y 2 

M^, My = Bending moments 

= Twisting or torsional moment 


(32) 
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Each of these forces influence the yield behavior at a cross section and the 
inelastic response depends upon the interaction between them. In the past, 
elliptical, parabolic, spherical (ref. 2) and other forms of yield functions 
have been used in inelastic analysis of structures. In the present study 
a spherical yield function is used as indicated below 


f* = 


(^) 2+ (^) 2+ (^) 2 + (£-) 

\ ox / ' ox / ' oy ' ' oz / 


= Y 


(33) 


in which Y denotes the normalized yield value which may change during straining 

and P , M , M and M are the characteristic values of axial force, 
ox ox oy oz 

torsional moment, and bending moments about y and z axes, respectively. 

Similar expressions can be written for a plate element (ref. 1). 


The above yield function can be used in conduction with an average force 
model (refs. 1, 3) which assumes that an element undergoes plastic deformations 
if the loading function f* determined from average values of stress resultants 
acting at the member ends exceeds the current normalized yield value. 


Flow Rule 


A flow rule expresses the relationship between plastic strains and 
stresses. In the present study rather than solving the flow equations 
rigorously in terms of stress resultants, an approximate procedure (refs. 1, 3) 
is employed. In this method if f* is the plastic potential at the end of 

a load increment in dynamic analysis (computed using member forces P. 

est 

obtained from an elastic analysis) , then the increment of plastic potential 
df*^ f° r an element already undergoing inelastic deformation is obtained from 
the equation 


df 


* 

1 



* 

f 

prev 


(34) 


* 

where fp rev is the plastic potential at the end of previous load or time 

increment. For a transitional element the equation for the increment of 
plastic potential is of the form 


df 


1 


* 


* * 
est o 


where f * is the initial yield value, 
o 


(35) 
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As shown in detail in reference 1, the correct value of the plastic 
potential increment df* can be computed from the relation 


* 

df - 


K, 


1 + K 


o * 

— df i 


(36) 


The current corrected value of the plastic potential f is then obtained as 

curr 


* * 
f “ f 4* df 

curr prev 


(37) 


for an element already undergoing plastic deformations and 

* * 

f = f + d * 
curr o f 


(38) 


for an element entering the plastic range. 

* 

With the known value of the plastic potential f , the final values of member 
(element) forces are computed by a proportioning procedure represented as 
(refs. 1, 3) 


?. = P. 

1 1 ^ 
curr est 


curr 


‘est 


(39) 


The overall numerical procedure utilized in the computation of the element 
and structural responses is the same as the one outlined in references 1 and 2. 


NUMERICAL RESULTS 


To determine the feasibility of the proposed method, two structures for 
which results are available in the literature were analyzed and the results 
were compared with the works of other investigators. 

Example 1 - The first example considered consists of the simply supported 
beam with the I cross-section and material property, as shown in figure 1. 

The given stress-strain curve corresponds to material B as defined by 
Chajes (ref. 4) who had originally studied this structure. As in the case of 
the above reference, the stress-strain curve is idealized by a bilinear 
relationship. This is then utilized to obtain the moment curvature and the 
normalized universal stress-strain curves, as depicted in figure 2. 
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Subsequently, the response of beam when subjected to either a concentrated 
load at its midspan or a uniformly distributed load over its entire length 
is obtained. It should be mentioned that Chajes (ref. 4) has presented 
a "closed-form" solution for the deformation characteristic of the midspan 
of the beam when subjected to a concentrated load. However, his solution is 
based on the assumption that only the flanges resist the bending moment. In 
figure 3 are shown the deflection responses of the midspan as obtained in this 
study as well as that reported in reference 4. As can be observed, close 
agreement between the two sets of data is exhibited. 

Example 2 - As a second example, a pin-based portal frame studied, both 
theoretically and experimentally, by Takahashi and Chiu (ref. 5) is analyzed. 
The structure consists of W12X27 sections, arranged to deform around their 
strong axes. The geometry of the structure and its loading are shown in 
figure 4a. The idealized structure and the equivalent nodal loading is shown 
in figure 4b. Note that the effect of the girder depth has been taken into 
account by introduction of an additional equivalent load and a bending 
moment at the top of the loaded column. In the analysis, a curve is fitted to 
the stress-strain relation for the mild steel with an average yield stress of 
255 MN/m^ and modulus of elasticity of 203.5X10^ MN/m^. This is then used 
as a basis for determination of moment-curvature relationship for the structure. 
In figure 4c is shown the horizontal deflection response of the top of the 
unloaded column. Again, close correlation is observed. 
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